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Abstract 

We study H— theorems associated with the Brownian motion with constant drift on the hyperbolic 
plane. Since this random process verifies a linear Fokker-Planck equation, it is easy to show that, up 
to a proper scaling, its Shannon entropy is increasing over time. As a consequence, its distribution 
is converging to a maximum Shannon entropy distribution which is also shown to be related to the 
non-extensive statistics. In a second part, relying on a theorem by Shiino, we extend this result 
to the case of Tsallis entropies: we show that under a variance-like constraint, the Tsallis entropy 
of the Brownian motion on the hyperbolic plane is increasing provided that the non-extensivity 
parameter of this entropy is properly chosen in terms of the drift of the Brownian motion. 



*Electronic address: vignat@u-psud.fr 



1 



I. INTRODUCTION 



The existence of an if— theorem for a statistical system submitted to some constraints 
ensures that the Shannon entropy h (ft) of its probability density function (p.d.f.) f t at time 

t, 



is increasing with time. As a consequence, the asymptotic (stationary) p.d.f. foo (x) of the 
system is the maximum Shannon entropy p.d.f. that satisfies the given constraints. 

The Shannon entropy is a particular member of a family of information measures called 
Tsallis entropies that were introduced Q in 1988 by Tsallis i the context of statistical physics; 
they are defined as 



where q > is the non-extensivity parameter. It can be checked using l'Hospital rule that 
the Shannon entropy is the limit case lim^i h q (ft) = h (ft) . 

A natural question then arises: under what conditions does an H— theorem extend to an 
^.-theorem, where the Sha„,to„e entropy is replaced by a Thallis entropy 7 Several studies 
have been devoted to this problem in the recent years: for example, Plastino et al [3] study 
the conditions of existence of an H q — theorem for the following non-linear Fokker-Planck 
equation 



for some parameter q, while Tsallis and Bukman [2] apply the same approach to the equation 



for some positive parameters [i and v. 

Our aim in this paper is to show the existence of both an H— and an H q — theorem 
for the Brownian motion - more precisely its x— component - with constant drift on the 
hyperbolic plane. This study is simplified by the fact that this component of the Brownian 
motion satisfies a linear Fokker-Planck equation, for which the conditions of existence of 
//—theorems are well-known, as we will see below. However, even in this simple case, this 
study reveals an interesting link between the constant drift parameter of the Brownian 
motion, the constant negative curvature of the hyperbolic plane and the non-extensivity 
parameter q that characterizes the entropy. 
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II. THE CLASSICAL ENTROPIC APPROACH TO THE LINEAR FOKKER- 
PLANCK EQUATION 



A. General approach 

In the case of a system described by a univariate p.d.f. f t (x) that verifies the linear 
Fokker-Planck equation 

= ~ (K (x) f t (x)) + -^(Q (x) f t (x)) , (III) 

where K fx) is the drift function and Q (x) the diffusion function, it can be shown (see for 
example 8]) that the relative entropy (or Kullback-Leibler divergence) 

h(f t \\g t )= [ f t (x)\og^\dx 
Jr 9t {x) 

between two any solutions of ( III.ip decreases to with time. More precisely, it holds 

|M A ||») = - 1 Q M /, W (A log £g) 2 < o (II.2) 

since the diffusion function is assumed positive. Thus entropy is decreasing and bounded, 
so that the limit distributions and verify 



i Q{x)ft{x) (^ l0g 



ft(x) x 2 



dx = 0, 



9t \ x ) 

which implies that f^ and g^ coincide. This proves the unicity of a stationary solution of 
the Fokker-Planck equation fill. II) . 



We note that the quantity 



is nothing but the relative Fisher information [ldl . eq. (174)] between the p.d.f.s / and g; 
thus the integral that appears on the right-hand side of (1II.2P is a weighted version of this 
relative Fisher information, with the diffusion function Q (x) as the weighting function. 

In order to deduce an H— theorem from this result, we denote as g^ the stationary solution 
to (III.ip . assuming that it exists. We then remark that the relative entropy between any 



solution ft of ( III.ip and the stationary solution g^ is related to the Shannon entropy of f t 



MS 



h(ft\\goo) = -h(ft) - / ft (z) log g^ (x) dx. 



Thus, provided that the solution f t verifies at any time the constraint 

/ f t (x) log (x) dx = r?i (II.4) 

JR 

where rji is a constant, we deduce that as the relative entropy decreases to 0, the Shannon 
entropy of the solution f t increases with time to its maximum value. 

We remark that, in a statistical physics framework, the relative entropy h{f t \\g 00 ) coin- 
cides with the free energy. Moreover, taking the limit t — > +00 in (III. 4ft shows that the 
constraint 771 is also equal to the negentropy of the stationary solution g^, namely 



/ 

Jr 



g^ (x) log goo 0) dx = -h (goo) • (H.5) 



B. The x— component of the Brownian motion in the Poincare half-upper plane 

In jlj], Comtet et al. derived the differential equation verified by the x— component X t of 
the Brownian motion with constant drift /j and constant diffusion constant D in the Poincare 
half-upper plane representation of the hyperbolic plane, 



This is a linear Fokker-Planck equation; in the notations of (IH.ip . the diffusion function 
is quadratic and positive, Q (x) = D (1 + x 2 ) and the drift function is linear, K (x) = 
D (1 — 2/i) x. Moreover, for a positive drift /i, the asymptotic solution reads 

9oo(x) = A il {l + x 2 Y^ (II.6) 

with a normalization constant A, = F ^ + / 1 \ . The constraint to be verified by the p.d.f. f t 

v r(/4)r(|) 

is thus deduced from the above results as 

/ f t (x) log (1 + x 2 ) dx = 771 Vt > 0. (II.7) 
Jr 

(note that the normalization constant log and the constant —fi — ~ that appear in log goo 
need not be taken into account). 

The value of the constant 771 can be easily computed as 



'/j - '' ( fi+ ^ ) - ' '('/') 



where ip is the digamma function. 



We note that this constraint can be imposed by a simple scaling of the random process 
X t since the function 

a h-> ft (x) log (l + ax 2 ) dx 
Jr 

is a bijection from IR + to IR + . We also remark that this scaling does not require the existence 
of a variance for X t ; for example, the Cauchy p.d.f. 

/c(x) = -— xeR 

7T 1 + X Z 

- which is the asymptotic distribution of the Brownian motion without drift (/i = 0) on the 
hyperbolic plane - has an infinite variance but finite nonlinear moment 



/ f c (x) log (1 + x 2 ) dx = 2 log 2. 
Jr 

We deduce the following theorem. 



Theorem. The Shannon entropy of the x— component, normalized according to ( 171. 7j) . of 
the Brownian motion on the Poincare half-upper plane representation of the hyperbolic plane 
increases over time. 

III. A GENERALIZATION TO THE TSALLIS ENTROPIES 
A. General approach 

The monotone behavior of the relative Shannon entropy between any two solutions f t 
and g t of the linear Fokker-Planck equation (IH.lj) has been extended by Shiino [3] to the 
case of the relative Tsallis entropy, defined as 

K (ft\\9t) = -^r ( [ f (x) g l - q (x) dx - 

1 ~ 1 \JR 



More precisely, the derivative with respect to time of this relative entropy verifies 11] 



^h q {f t \\g t ) = -q [ D(x)f t (x)(^) 9 (|-bg£M) dx < 0, Vg > 0. (TTT.l) 
dt J R \9t{x)J \ox g t (x 



We note that this inequality holds for any positive value of q and simplifies to (III. 2ft as 



q — > 1. Moreover, the right-hand side integral in ( 1III.1I) can be considered as a q— version of 



the relative Fisher information that appears in f ill. 21) 
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In order to deduce from this monotonicity an H q — theorem, we need the additional as- 
sumption that the stationary solution g^ can be written under the form 

2oc (x) = C q (l + U(x))^ 

for some specific value q = g* of the non-extensivity parameter. This assumption means 
that the p.d.f. is itself a maximum g* — entropy pdf with constraint 

J g%(x)U(x)dx = r lq , {111.2) 



for some constant value rj qt . We refer the reader to [6] for the conditions on the diffusion and 
drift functions of the Fokker-Planck equation that ensure the validity of this assumption. In 
this case, we have 

K (/tltooo) = -<%*K (ft) - t~~ f ft (*) u (*) dx ~ P* 

where /3g„ = — ^ is a constant. Thus, provided that the constraint 

ft (x) U (x) dx = r] qt 

is met at all times, we deduce that the g*— entropy is increasing with time and reaches 
asymptotically the entropy of the stationary pdf g^ . 



B. The x— component of the Brownian motion in the Poincare half-upper plane 

In the special case of the x— component X t of the Brownian motion in the Poincare half- 
upper plane, the stationary solution flll.6j) verifies the condition fjIII.2j) with the value g* 
such that —fi — | = j^— and the function U (x) = x 2 , namely 



1-g 



x ft (x) dx = r) qt (III. 3) 



with 



Q* = 1 < g* < 3. (III.4) 

As a consequence, since the g— relative entropy is decreasing, the q— entropy is increasing 
under the conditions that q is chosen equal to g* as in (lIII.4j) and that the constraint (1IH.3|) 
is verified. 
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Moreover, by taking the limit as f -> +00 in (1III.3j) . the constraint rj qt in 3 j) is 
computed as 

% = ~, ~ / x 2 g% (x) dx. 

The value of this constraint is 

g* 



1 V ( 



2 



We note that this constraint can always be imposed by scaling: denoting ol (X t ) = 
(1 — q)^ 1 J x 1 fl (x) dx the value of this constraint and cr 2 {Y) the one corresponding to 
Y = aX with a > 0, then 

a\ (Y) = a 3 -V| (X) . 

We deduce the following result. 



Theorem. With q = q* as in (III. 4}) , the q— entropy of the Brownian motion, normalized 



according to UII 3\) . is increasing with time. 



IV. NUMERICAL ILLUSTRATION 



The following figures depict ten realizations of a discretized version of this Brownian 
motion, with parameters D = 0.01, m = and \i = 1 on Figure [1] while D = 0.01, m = 3 
and fx = 7 on Figure. In both cases, the process starts from the point x — 0, y — 1 and 
the superimposed thick curve is the asymptotic probability density of its x— component. 
Without external drift (Figured]), the random process wanders for a long time far away 
from the real axis before "falling" on it. The distribution of the "landing points" on the 
real axis is thus very wide, in fact a Lorentz distribution with infinite variance. With an 
external drift (Figure [2]), the process is forced to walk in the direction of the real axis, so 
that the landing points are more concentrated around 0, what is reflected by their narrow 
distribution, a q— Gaussian distribution with variance a 2 = 0.2. 
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Figure 1: ten realizations of a discretized version of the Brownian motion on the hyperbolic plane 
with parameters D = 0.01, m = and fi = 1 

3 -°i 

2.5- 
2.0- 




Figure 2: ten realizations of a discretized version of the Brownian motion on the hyperbolic plane 
with parameters D = 0.01, m = 3 and /i = 7 

V. THE BROWNIAN MOTION IN THE UNIT DISK 

Another representation of the hyperbolic space is the unit disk © 
= re %e , r < 1, < 9 < 27r} with the metric in polar coordinates 

ds 2 = ^ (dr 2 + r 2 d9 2 ) . 

There is a conformal mapping between the Poincare upper half-plane H 
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{z = x + iy, y > 0} and the unit disk D denned as 

iz + 1 
w = . 

Z + I 

Comtet et al show that the density of the radial component of the Brownian motion in the 
unit disk representation of the hyperbolic space converges to 6 (r — 1) as t -> +oo. Having 
no explicit Fokker-Planck for the radial part 9 t of this process, we were unable to prove a 
corresponding H— theorem. However, we show here that we can use the maximum entropy 
approach to derive the asymptotic distribution of the angular part 9 t , using the following 
result: 

Theorem. If the random variable X has maximum entropy under the log- constraint 
E\og (1 + X 2 ) = 7, then the random variable 

x x 



has maximum entropy under the constraint E\og (l — X 2 ^ = —7. 
More precisely, if the p.d.f. of X reads 

then E\og (1 + X 2 ) = ift (fi + |) — ^ (/f) and the p.d.f. of X reads 

f ( £) = Lfr±il (1 _ zy- 1 

JXK ! r(|)r(/i) v ; ~ - 

with E\og (l - X 2 ) = -^log (1 + X 2 ) = V (/i) - ij) (n + \) ■ 

Since in the asymptotic regime, f r (r) = 5 (r — 1) , the conformal mapping becomes 

2X . X 2 -l 

cos6 = TTx- 2 > S1 ^ = ^TT' 

it can be easily verified that 

(0 n\ X ~ 

cos = . = X 

\2 Aj VTTx 1 

so that a simple change of variable yields 

6~ (l-sinflf- 1 

as obtained in 

The asymptotic distribution of the angular part has thus maximum Tsallis entropy with 
parameter q such that 

q = - — 7 < 1. 

(i- 1 
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VI. CONCLUSION 



We have shown how to use the monotonicity of the Shannon or Tsallis relative entropies 
to deduce an //—theorem for the Shannon and Tsallis entropy, first in the general case 
and then in the case of the Brownian motion on the hyperbolic plane. Three remarkable 
results have been observed in this study: first, the natural drift induced by the negative con- 
stant curvature of the hyperbolic plane transforms the asymptotically Gaussian of the usual 
Brownian motion on the plane to the Cauchy distribution, which belongs to the extended 
family of Tsallis distributions. Secondly, the addition of a constant positive external drift 
transforms this Cauchy behavior to a non-extensive behavior with non-extensivity parame- 
ter directly related to the value of this drift. This appearance of non-extensive distributions 
in the context of an underlying curved space remains to be linked to physically relevant 
experiments and data. At last, the Tsallis distributions with q > 1 on the Poincare half 
upper-plane realization of the hyperbolic plane transform, via the conformal mapping, into 
Tsallis distributions with q < 1 on the unit disk realization of the hyperbolic plane. 

VII. ANNEX: PROOF OF SHIINO'S RESULT 

Following Shiino's notations, we consider 

D q (ft\\9t)= [ ftixygtixf-^dx 
Jr 

and omit the time and space variables for readability. The time derivative reads 

|d. (fh) = 9 / g/) + a- 9 )/ /v« (| 9 ) . 



Since / and g are both solutions of the Fokker-Planck equation flll.ip . we deduce 

1-1 / Q Q2 

dx dx 2 



d 2 

(Qg) 



dx 2 

The integrals with the drift function K (x) can be integrated by parts, yielding respectively 
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and 

so that their sum vanishes. 

The integrals with the diffusion function Q (x) are also integrated by parts according 
respectively to 

f\ rl d 2 r d f fV^ 1 d 



g J dx 2 J dx \g J dx 



and 



g J dx \ g J dx 



g J dx \g J dx 

~ (l-q)q fQ(-) q l -(-) — 
J \g) dx \gj dx 

Their sum consists in an integral with the diffusion function Q 

^ ( g ) ( dx ( g ) da: da; ( ^ g ) da; 



(q-l)qj 
and an integral with its derivative 



.<?/ \ da; \ (?y dx \g J J dx 
The second integral is easily seen to vanish, while the first one can be simplified to 
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Thus the Tsallis divergence h q (f\\g) = ~zp[D q (f\\g) verifies the stated equality. 
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